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The dynamic model construction of transmission network components that 
include generator buses, load buses and power branches, within MATLAB- 
Simulink environment is presented. The degree of frequency deviation of 
buses when the power of motor loads and static loads vary, is studied. 
Furthermore, the influence of motor loads with different inertia constants are 
considered in the control technique of load frequency using a PID controller. 
The results show that the frequency oscillation of the dynamic load is greater 
than the frequency oscillation of static load. Also the speed of frequency 
control of the dynamic load is greater than the speed of the frequency control 
of the static load and the inertia constants of the dynamic load has significant 
influence on the frequency control characteristics. 
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1. INTRODUCTION 

In recent years, the increasing stress in a transmission system may limit the effective power transfer 
from generation to the load. At steady state, the frequencies on different buses in a system are synchronized 
to a common nominal value and the mechanical power is balanced with the electrical power on each bus. 
Suppose a small change in power balance injection occurs on an arbitrary bus. This mismatch between 
generation and load causes the bus frequency to deviate from its nominal value. Consequently, the system 
may lose stability and cause damage to the facilities if the frequency deviations are not tightly controlled 
around zero. This is important as the frequency of a power system is an indispensable performance signal to 
the system operator for security and stability considerations. Hence, the desired power system frequency 
should be kept within a very small, acceptable margin around its nominal value. Otherwise, the system 
operator must take relevant and effective actions immediately to avoid physical damage to devices. 

A power system is a combination of generation, transmission, distribution networks and loads. The 
active and reactive power demands from different loads vary continuously. The induction motors account for 
a large portion of electrical loads in industry and air conditioning in business and residential areas. Dynamic 
characteristics of these loads affect the response of system stability when compared to static loads. Hence, the 
influence of motors as dynamic loads on transient stability studies have attracted a great deal of attention. 
The linearized models of dynamic loads [1-3] and load frequency control strategies [4] on a nonlinear system 
are considered to develop a robust controller. 

The real dynamic loads are composed of a variety of motors whose parameters vary over a large 
range, so the aggregation of motor load with suitable parameters, especially inertia constant H is very 
important. This constant significantly influences the transient and dynamic system behavior [5]. In this 
context, the presence of motor load in a network with influence of inertia constant to frequency response is 
investigated and reported. 
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In an interconnected system, as a power load demand varies randomly, both bus frequency and line 
power interchange will also vary. The objectives of load frequency control (LFC) are to minimize the 
transient deviations in these variables (bus frequency and line power interchange) and to ensure their steady 
state error to be close to zero. Thus, a control system is essential to overcome the effects of the random load 
changes and to keep the bus frequency within the range of nominal value. 

In the past decades, many researchers have investigated the power network frequency regulation. 
Since the system frequency is essentially related to real power balance, it is natural to control real power 
output in the generation side through the use of automatic generation control (AGC) system [6-9]. 

There are some simulation studies of frequency load control in a single area of power network 
which include one generator or several generators in parallel operation [10-11]. However, for large scale 
power systems which consist of interconnected control areas, it is important to keep the frequency constant 
and the tie line power near to its scheduled value. This is achieved by implementing LFC as Tie line 
controller and PID controller [7-8]. There are much more analytic studies that relate the behavior of the loads 
and frequency responses of interconnected power network [12-14]. This is normally achieved through 
establishment of the well-known generation swing equation and implementation of classic tools such as 
transfer function [12-13] and state space method [15] in the swing equation to estimate the amount of 
frequency deviation. 

With the above motivations, this paper contributes initially to the construction of a dynamic model 
of a transmission network including generator buses, load buses and power branches within MATLAB- 
Simulink environment. Then investigates the degree of frequency deviation of buses when power of motor 
loads and static loads vary. In addition, the influence of motor loads with different inertia constants are 
considered in the control technique of load frequency using a PID controller. 

The remaining parts of this paper are organized as follows: mathematical modelling of a dynamic 
system are shown in Section 2; the proposed controller of load frequency and influence of motor load to 
control quality are presented in Section 3; the simulation results are given in Section 4; and the concluding 
remarks are given in Section 5. 


2. NETWORK MODEL DEVELOPMENT 

In the interconnected network, buses could possess generators, loads or neither. A generator bus not 
only has an AC generator that converts mechanical power into electrical power through a rotating prime 
mover, but some loads may also be attached to it. However, on the load bus, only loads are connected to the 
load bus. 

To construct the dynamic model of a network, the description model of a Turbine-Generation 
system needs to be developed. Then the models of power equilibrium equation for the generation bus, load 
bus and power branch between buses need to be established. All of these are essential in the interconnection 
network for studying dynamics network. 

To proceed with these investigations, the following assumptions are considered in this paper: The 
lines are lossless and characterized by their reactance Xy. All voltage magnitudes of buses are equal to their 
nominal values. Reactive power injection on the buses and reactive power flow on the lines are ignored. 
Rating power of generators is similar in the network model. 


2.1. Mathematical modeling of turbine-generator system 

In terms of AGC under the change of load frequency, the conditions needed to be known are the 
network model, analysis of load characteristics, linearization of swing dynamics on generator buses, power 
flow dynamics on the branches, and a measure of dis-utility to users when they participate in frequency 
control. 

In an interconnected power system, AGC equipment is installed for each generator. In steady state, 
the change of mechanical power AP Mi will be a constant. Any change in load is reflected in the frequency. 

A study of the system for small changes around a nominal setting, the Turbine-Generator system 
may be represented by the time constant of turbine T t and the time constant of governor TG [9-15]. 

The generator response is considered to be instantaneous: 


dAP M j 

dt 


T ^ AP *Mi F kXgvi 


(i) 


dAXgvi 

dt 


T G i 


-Ax n 


T Gi R i T Gi 


( 2 ) 


Laplace transformation of (1) and (2) yield: 
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APmi( 5 ) — 


(l+sT t 0(l+sT G 0 


[AP ci (s)-^AF i (s)] 


(3) 


2.2. Mathematical modeling of power equilibrium equation 

The increment in power input to the generator-load system is defined as (A P Gi — A P Li ). 


AP Gi ~ APMi ^Ploss,Gi 


AP Gi ~ APMi ^Ploss,Gi 


(4) 

(5) 


This increment in power input to the system is accounted for in three ways [14]: 

(i) Rate of increase of stored kinetic energy in the generator rotor. At scheduled frequency (f°), the stored 
energy is: 


Kel = H ci X P ri 


( 6 ) 


The kinetic energy being proportional to square of speed (or frequency), at a frequency of (f°+Afi) 
the kinetic energy is given by: 

W kei = wg el 

Therefore, the rate of change of kinetic energy is: 


( /°+Af,-\ 2 

TT _ (. 2Af,\ 


( f° ) 

^Gi^ri (^4 + ~fo~j 

(7) 


dt 


(w kel ) 


f° dt 


(Aft) 


(8) 


(ii) As the frequency changes, the motor load changes, being sensitive to the speed, the rate of change load 
with respect to frequency: 


D Li Afi= ^Af t 


(9) 


where D^Af* is the load change that is frequency sensitive. D L i=-^jr (MW/Hz) expressed as % 

change in load divided by % change in frequency. D Li can be determined empirically. 

(iii) Rate of increase export of power lines AP Une i 

The power equilibrium equation could be expressed as: 

A Pci ~ AP U = (Af t ) + D Li Af £ + A P Unejl (10) 

Dividing (10) by P r and rearranging it in per unit (pu): 

hP G i( pu) - AP Li (pu) = Aft(pu)] + f>Li Af;(pu) + A P lineil (pu) (H) 

From (4) and (5) the change in the generator power is derived as A P Gi = A P Mi — D Gi Af t , which 
could be substituted in (11) to yield (12) and then (13) in pu: 

7 T^[Af;(pu)] + D Gi Af ; (pu) + £> Li Af;(pu) = AP Mi (pu) —AP line i (pu) - AP Li (pu) 

^fA f(pu)+ (D Ll + D G t)Af(pu) = AP Ml (pu)-AP ltne>l (pu) - AP Ll (pu) 

2.3. Mathematical modeling of generator bus 

Each generator bus could be represented as Figure 1 [16] where: 


( 12 ) 

(13) 


A*Wf(pu) = Apr(pu) - a P} n (pu) 


(14) 
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Figure 1. Power diagram of generator bus 


The corresponding power equilibrium equation in per unit could be expressed as: 
fi(pu) + (D Ll + D Gi )Afi(pu) = A P Mi (pu) ~ APf ut (pu) + AP, m (pu) - AP Li (pu) ( 15 ) 

Let D t = (D ci + D Li ); K pi = 1/D h T pi = ^ 

Laplace Transform (LT) of equation (15) yields: 

AF ^ s ) = (i - Ap ^ - Ap ° ut <»+ Ap / n ( s )] (16) 


2.4. Mathematical modeling of load bus 

The load on a power system consists of variety of electrical drives. There are either the dynamic 
loads (e.g. induction motor) or the static loads A dynamic load model that can expressed as in [10]: 


pf _ pf® 

r Loadi r LoacLi 


= M'Lai = D Li A ft + h(A.fo = 0 


Where D Li is the damping of load at bus i th 

The rotating masses of Induction Motor (IM) load following kinetic energy: 


(17) 


w kei (f) = \jt<of 


(18) 


Where HpJiCOi with Ji is moment of inertia 

The change in the kinetic energy, which is equal to the power consumed P^i by IM is given by: 

n _ dw kei (J) (19) 

Flm - dt 

And: 


_a p _ dAWfcg^(/) _ 2Hj]y[iP r i d /. r>, 

IMi ~ dt ~ f° dt ^ ) 

From (17) and (20), the change of motor power in pu as: 

^aat(P^=D L t A/i+^f^T 

A load at bus i th is modeled by the equation in pu: 


2 HlMi 


^A/iCpu) + D Li Afi(pu) = —AP L ,(pu) - AP° u£ (pu) + AP'”(pu) 


Let V • = — • K 1 ■ = — 

pl D Li /o’ K P' D Lj 

Laplace transform (LT) for (22) yields: 
K/ n 


AFi(s) 


(1+ T’ pi s ) 


[-AP Li (s) - AP° ut (s) + APt n (s)] 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 
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A static load model that can expressed as in [17]: 


A P, 


static _ 


Loadi 


D LiA/j 


(24) 


The frequency deviation A/* is calculated similar to induction motor but without change in store kinetic 
energy. 


D Li Afi(pu) = —AP L j(pu) - A/T'CPU) + AP/ n (pu) 


(25) 


2.5. Mathematical modeling of of branch power 

Assume that the line admittance parameters are purely imaginary, branch power connected to bus i t] 
is yielded: 


Pi =^|K i ||K fc |B ifc sin(0 i -0 fc ) 


(26) 


k=1 


Where: V t = \Vi\a6i , V k = \V k \/-6 k . In this case I V t \ \V k \B ik are constants. 

Linearizing Pi around the operating point: P t = P® + A P t and 9 t = 0° + A 6 t , 9 k = 0 k + A 6 k 

n 

Pi = P? + A P t = JV ; | Wk\ B ik sin(0° + - 9° k - A9 k ) 

k=1 
n 

= ^Vil Wk\Bik [sin(0° - 0j“) cos(A6»i - A 9 k ) + sin(A0i - A9 k )cos(9f - 9 £) 


(27) 


fc = l 


Let assume that the angle of increments tends to zero: 
cos(A Q t — A0 k ) « 1; sin(A0j — A0 fc ) « (A0j — A0 k ) 

n 

A= JV ( | |F fc |^ fc cos(0P - 0fc)(A0j - A0 k ) 

fc = l 

Let T' ik = \Vi\\V k \B ik cos(9i — 9 k ) 

n 

AP i = Y J T'ik{.M i -A9 k ) 


(28) 


(29) 


fc = l 


Assume coherency between the internal and terminal voltage phase angles of each generator so that these 
angles tend to “swing” together: A 6 t = A S t 


AP i (pu)=Y j T’ ik (AS i -A8 k ) 


(30) 


k=1 


n 

APiipu ) = 2n ^ T' ik (J A ftdt - J Af k dt} 

k=l 


(31) 


Let T ik = 2nT' lk . 

Laplace Transform (LT) for equation (29) yields: 

1 n 

AP;0 } = t ^ af ^ ~ ( 32 ) 

k=1 
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2.6. Dynamic network model 

In summary the dynamic model of the transmission network is specified by (15), (22) (25) and (31). 
To simplify notation, let drop the (pu) from the variables denoting per unit and write (15), (22) (25) 
and (31) as: 


yA/i + (D Li + D Gl Wi = AP M ; - AP° ut + AP[ n - A P Li 

(33) 

+ D Li Afi = -A P Li - A P° ut + A Pl n 

(34) 

D Li Afi = —AP Li - AP° ut + Aif 

(35) 

71 

A Pi = 2n ^ T' ik ( J Afcdt - J Af k dt ) 

(36) 


k=1 


The model (33)-(36) captures the power system behavior at the timescale of seconds. 


3. THE PROPOSED LOAD FREQUENCY CONTROLLER WITH CONSIDERATION OF THE 
MOTOR LOAD INFLUENCE 

In this paper, consider all loads at buses to be induction motors. They are compared to static loads 
used in the same network model. These dynamic loads will bring significant impact to the power system in 
term of frequency due to the kinetic energy that is stored in the rotating masses of the motors [18]. Dynamic 
load modeling mentioned in previous part is widely used as a model to analyse power system stability 
problems. The effect of the induction motor on the stability of the system is analyzed by changing inertia 
constant through small disturbance. A large inertia constant of motor load entails slow frequency dynamics in 
response to generation-load imbalance. 

As there is a step change in load, the turbine-governor takes action that changes the mechanical 
power injection in response to frequency deviation to rebalance power. It is possible to include any secondary 
frequency control mechanism such as AGC that operates at a slower timescale to restore the nominal 
frequency. A proposed two-area control scheme for LFC [17, 19] is applied to control frequency in this 
paper. The simulated responses of LFC model will be compared in the tie line controller (without PID block) 
[11] and the proposed PID controller given in Figure 2. This comparison is based on the oscillation and 
response time of frequency at buses in the network. This indicates the relative stability of the system. 

The frequency control model of the two area system but with the additional PID controller is shown 
in Figure 2. Suppose there is a change in load (AP Li ) in any bus in the system, the frequencies at buses in the 
system settle to a steady state value. 



Figure 2. Propose control block diagram of two-area interconnected system 
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4. RESULTS AND DISCUSSIONS 

4.1. Sub frequency deviation response in dynamic network model 

In this paper, the proposed approach is applied on a 230 kV power system of Figure 3. The system 
consists of three thermal generating units, the connected branches that have the same impedances, and the 
loads that are placed on buses 2, 3, 4 and 5. The frequency deviation response is illustrated with a sudden 
load increase of 10% (AP=0.1 p.u) on bus 3 at t = 0.2s. For the purpose of analysis, performance of the 
dynamic system model and static model under various scenarios are compared. The system data are given in 
the Table 1,2,3 and the Induction motors data are r| (efficiency) =80%; pf (power factor)= 0.85 lag; Pout = 
50HP; and inertia constant H IMi =6. It is assumed that static loads connected to the four buses (bus 2, bus 3, 
bus 4 and bus 5) have equal power magnitude as induction motors at each respective bus. 



Figure 3. Test system diagram 


Table 1. Generator Data 


No.Gen 

Type 

H 

R 

X’ d 

T g 

T t 

D g 

1 

Steam 

9 

0.03 

1.3125 

0.2 

0.3 

0.1 

2 

Steam 

9 

0.03 

1.3125 

0.2 

0.3 

0.1 

3 

Steam 

9 

0.03 

1.3125 

0.2 

0.3 

0.1 


Table 2. PID Controller Data 


No. PID 

K P 

Ki 

K d 

1 

15 

10 

7 

2 

20 

15 

10 


Table 3. Bus Data 


Bus 

1 Vi 1 (pu) 

5i 1 (pu) 

D l 

1 

1.017 

-1.02 


2 

1.0566 

-2 

0.6 

3 

0.996 

-4 

0.6 

4 

1.013 

-3.7 

0.6 

5 

1.0566 

-2.2717 

0.6 


The test network is modeled in the MATLAB/Simulink environment. Several scenarios of 
connections among generator buses, load buses and power branches listed in Table 4 are considered. The 
loads are modelled as frequency dependent loads. The responses of the frequency deviations at bus 1 
(generator 1) are chosen as a reference bus for comparison between motor loads and static loads under 
various scenarios as shown in Figure 4 to Figure 13. 
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Table 4. Cases of Changing Connected Networks 


Case 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

G1 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

G2 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

G3 







J 

J 

J 

J 

Line 1-2 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

Line 1-3 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

Line 2-4 


J 

J 

J 

J 

J 

J 

J 

J 

J 

Line 2-3 



J 

J 


J 

J 

J 


J 

Line 3-4 




J 

J 

J 

J 

J 

J 

J 

Line 3-5 






J 


J 

J 

J 

Line 4-5 







J 


J 

J 



Figure 4. Frequency deviation in case 1 Figure 5. Frequency deviation in case 2 



Figure 6. Frequency deviation in case 3 


Figure 7. Frequency deviation in case 4 




Figure 9. Frequency deviation in case 6 
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Figure 12. Frequency deviation in case 9 



Figure 13. Frequency deviation in case 10 


Increasing the load by 10% at bus 3, in the case 1 there is small drop in frequency deviation at bus 1 
and steady state response time is almost 27 seconds in both static loads and motor loads. However, the 
frequency deviation at bus 1 increases graduatelly for cases 2 to 10. When adding line 2-4 to network model 
in the case 2, there is a significant change of frequency at bus 1 and it takes more than 30 seconds to reach 
the steady state. In the case 3, the oscillation of frequency at bus 1 reduces compared to the case 2 due to 
there are one more line 2-3 connected to network. In addition, the frequency osscillation and time response 
are higher for cases 4 to 10. These changes are related to addition of more loads and lines at buses in the 
network. In general, more incremental oscillations of the frequency deviation response observed at buses 
with motor loads than the static loads. 

4.2. Influence of motor load to the proposed load frequency control response 

The proposed load frequency control system with the addition of a PID controller is implemented at 
generators G1 and G2 in the full network (case 10) where the total load at bus 3 is increased by 10%. The 
network parameters and the gains of the PID controller are given in the Appendix. In order to evaluate 
reliability of the control system, the frequency response of the proposed controller is compared with the 
conventional Tie line controller for bus 1, bus 2 and bus 5 ( generator buses) as shown in Figure 14-16 
respectively. Based on these results, the PID controller give the frequency responses less oscillation and more 
precise than Tie line conntroller. Therefore, the load frequency control system with PID controller is 
continued to implement comparasion freqency responses in the system between motor loads and static loads. 

The analysis of Figure 17-19 show that the speed of frequency control response of motor loads is 
faster than the frequency control response of static loads. This is due to among the motor loads and 
generators transfer kinertic energy together, they support together and help for frequency responses reaching 
the steady state quickly. 
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Figure 14. Frequency response at bus 1 



Figure 16. Frequency response at bus 5 



Time (sec) 

Figure 18. Frequency response at bus 2 



Figure 15. Frequency response at bus 2 



Figure 17. Frequency response at bus 1 



Time (sec) 

Figure 19. Frequency response at bus 5 


Finally, the different load inertia constants =(5; 5.5; 6; 6.5; 7) at bus 3 (increase 10% load) are 
implemented in the proposed PID control system to gain frequency responses at respective buses in the 
network. The simulated results are shown in Figure 20-22. The tentative analysis of the obtained results 
suggest that value of inertia HjMi of the motor has a direct effect on the initial slope and the time of the peak 
response of frequency. However, HjMi does not affect the final steady state value of frequency. 
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Time (sec) 


Figure 20. Frequency response in bus 1 


Time (sec) 

Figure 21. Frequency response in bus 2 




Time (sec) 


Figure 22. Frequency response in bus 5 


5. CONCLUSION 

The model of a dynamic network for study into the effect of a small change in load on deviation of local 
frequencies of the network is presented. The network is composed of generators, load buses and branches 
between the buses. The proposed dynamic model of the network for analyzing the frequency deviation 
response is implemented within MATLAB-Simulink environment to compare the effect of static loads and 
induction motor loads. The proposed dynamic network model can be developed in an interconnected network 
with multi generator buses, load buses and connected branches. This paper also demonstrates the reliable 
operation of the PID controller for frequency response and the influence of an induction motor through 
inertia constants on the characteristics of the proposed load frequency controller. The drawbacks of these 
PID controllers are fixed control gains while the system changes continuously in time. Also in the frequency 
stability study, the effects of load dynamics and especially induction motor loads on frequency recovery 
phenomena need to be considered more seriously. 
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